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-«^ I Abstract. It is now well known that a primordial chemistry, involving light elements produced during the 

nucleosynthesis period, might develop during the hydrogen post-recombination era. In particular, molecular 
ions and primordial molecules such as H2, HD and LiH will be produced. We summarize this primordial 

r^ ' chemistry after the recombination epoch, and then present a simple gravitational collapse model of a cloud. 

D I The potentiality of fragmentation of this collapsing protoclouds through the thermal instability is also discussed. 

|JL^ , We suggest that this study could also be extended to the CO molecule, because the carbon reservoir molecule 

_^ ■ CO has already been observed in high redshifts objects. 

I INTRODUCTION 

•Q ■ The study of primordial chemistry of molecules addresses a number of interesting questions pertaining to the 

^~; I thermal balance of collapsing molecular protoclouds. In numerous astrophysical cases molecular cooling and heating 
j^ ■ influence dynamical evolution of the medium, see for example the review done by Dalgarno and Mc Cray (1972) 
Qs I who present molecular cooling and heating in HI regions, or Shu (1997) who discusses the influence of molecules in 
Q^ i star formation. 

The existence of a significant abundance of molecules can be crucial on the dynamical evolution of collapsing objects. 
Because the cloud temperature increases with contraction, a cooling mechanism can be important for the structure 
formation, by lowering pressure opposing gravity, i.e. by allowing continued collapse of Jeans unstable protoclouds. 
O I This is particularly true for the first generation of objects. Many authors have examined this problem (Bodenheimer 
1968, Hutchins 1976, Palla et al 1983, Villere & Bodenheimer 1987, Tegmark et al 1997) introducing molecular 
coolants, mainly H2 molecules. From Fall and Rees (1985) we know that H2 molecules are important cooling agents 
and can lead, sometimes, to the fragmentation process of the object and thus provide a possible scenario of globular 
cluster formation for collapsing galaxies. 

Recently Puy and Signore (1997), Bougleux and Galli (1997) showed that during the early stages of gravitational 
^ ' collapse, for some collapsing masses, HD molecules were the main cooling agent when the collapsing protostructure 
5^ I had a temperature of about 200 Kelvins. In our last papers (Puy and Signore 1997, 1998) we have numerically 
calculated the cooling functions for H2, HD and LiH molecules by considering the excitation of the twenty first 
rotational levels. H2 is an homonuclear molecule, so radiative transitions within the lowest electronic state involve 
only quadrupole transitions, A J = 2. Between two rotational levels, A J = J + 2 — J the energy of the transition is 
given by 

Ej+2,j = 2kBr{2J + 3) = kTj+2j (1) 

where Tj^2,.i is the temperature of the quadrupolar transition J + 2 -^ J, Br the rotational constant (for H2, 
Br = 85.33 Kelvins) and k the Boltzmann constant. Due to the non-zero dipole moment the transitions involve 
dipolar transitions, A J = 1. Thus the energy of the transition J ^ J + 1 is given by: 

Ej+ij = 2kBr{J + l) = kTj+ij (2) 

Tj+i^j is the temperature of the dipolar transition J + 1 -^ J, Br = 64.3 Kelvins for HD and Br = 10.7 Kelvins 
for LiH. 

In this paper we consider a collapsing molecular protocloud made up only of H2 and HD molecules. We do not 
consider collisions with H atoms because we consider only the central part of the cloud where the molecules are 
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present. In this picture H is less heavier and is confined in the circumcloud. Moreover, we do not consider the 
evolution of abundances. Puy & Signore (1998) showed that the evolution of H2 and HD abundance was low, and 
abundance of LiH was negligible (sec also Stancil et al. 1996). Thus, the constant ratio 77 = nH2/nHD (where nH2 
and uhd are respectively the density of H2 and the density of HD) around the primordial ratio 
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nH2/nHD ^ 1430 



(in the framework of the Standard Big Bang nucleosynthesis, see Puy et al 1993). 

In section 2, we recall the primordial chemistry. In section 3, with these approximations, we analytically calculate 
the population of the ground state level and first state level in order to evaluate the molecular cooling and molecular 
heating due to HD and H2 molecules. In section 4, we estimate the molecular cooling during the collapse of 
protoclouds. Then, in section 5, the possibility of thermal instability during the early phase of gravitational collapse 
is discussed. Finally, in section 6, we suggest that this study could be extended to CO molecules which have been 
recently detected at redshifts z > 4 (Ohta et al 1996, Omont et al 1996, Guilloteau et al 1997). 

II PRIMORDIAL CHEMISTRY 

In standard Big Bang theory, chemistry took place around the epoch of recombination. At z^ 1000 the chemical 
species essentially were (see Signore & Puy 1999): 

H,H+,D,D+,He,Li 

As the Universe expanded and cooled, different routes led to molecular formation. Many authors have described 
in detail primordial molecular formation and evolution; in particular, Puy et al. (1993) have simultaneously solved 
the coupled chemical equations and the density and temperature evolution equations. Finally with initial conditions 
given by the Standard Big Bang Nucleosynthesis (Signore & Puy 1999), we obtain the following abundances (see 
Puy & Signore 1999): 

H2/H - 8.4 X 10"'^ 



HD/H - 10 
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LiH/H - 10"^^ 
In the following due to the very low abundance of LiH we will consider only the primordial molecules H2 and HD. 

Ill ANALYTIC CALCULATIONS OF COOLING AND HEATING FUNCTIONS 
OF H2 AND HD MOLECULES 

We consider here only the transition between the ground state and the first state. Due to the existence of a 
permanent dipolar moment for HD {jjlhd = 8.310"'* Debye, Abgrall et al 1982), the transition concerning HD 
molecule is a dipolar transition whereas for H2 molecule the transition is quadrupolar (because it has no permanent 
dipolar moment). 

A HD molecule 



The excitation of the ground state is radiative or coUisional, the equations of evolution concerning the populations 
{X^^ for the ground state J = 0, X^^ for the first excited state J = 1) are given by: 
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where Cq^ and C^^ arc the cohisional coefBcients characterizing the colhsions between HD and H2. In the 
approximation of MaxwcU-Boltzmann distribution wc have: 
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with T^(f ~ T^^ the transition temperature between J = -^ J = 1 , Tm the temperature of the matter inside 
the cohapsing cloud, B^^ and B^^ are the Einstein coefficients. The radiative density is 
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where i^iq and Tf^ are the corresponding frequency and temperature, h the Planck constant, c the velocity of the 
light and Tr the temperature of the background radiation. 
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characterizes the Einstein coefficient for the spontaneous emission (see Kutner 1984) and Bhd = 64.3 Kelvins the 

rotational constant (see Herzberg 1950). 

We have the following normalisation for the populations: 

Xf ^ + Xf ^ = 1 



Moreover we consider a quasi-instantaneous transition so 

dt ^ dt 
Thus from these approximations we find for the populations X^^ and X^^: 
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and 
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The probability of coUisional de-excitation is given by 
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and the probability of radiative de-excitation: 
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Thus, we can calculate the molecular cooling (collisional excitation followed by a radiative de-excitation): 

Ahd = nHDXrC^.FPr"''E^f (9) 

where E^^ = kT(^(P is the energy for the transition J = ^ J = 1. 

The molecular heating (radiative excitation followed by collisional de-excitation) is : 

r"^ = UHDX^B^fu^fPrE^^o'' (10) 

If we take into account the following relation between the Einstein coefficients: 
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then we deduce for the molecular cooling 
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and for the molecular heating 
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Notice that calculating the ratio molecular cooling to molecular heating we obtain the following simple expression 
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This relation gives us an important information concerning the thermal influence of the thermal function '^^^ = 
pff-D — A^^ during a gravitational collapse. When the matter temperature Tm is higher than the radiative tem- 
perature, we have i]^^ > 1, i.e. A^^ > r^^ so molecules cool the collapsing cloud. Thus during a gravitational 
collapse, the matter temperature of the collapsing cloud is higher than the radiative temperature of the cosmological 
background (let us recall that in our case we do not consider some external sources like quasars, first stars etc.). 
Thus in this study, we consider only the molecular cooling which is the dominant term in the cooling function as we 
have shown for a collapsing cloud (Puy & Signore 1996). Moreover if we compare the temperature of the dipolar 
transition for HD, J = -^ J ~ I, T^^ =- 128.6 K with the temperature of the cosmological radiation between 
the redshifts [5,20], i.e. T^ G [13.5,54], we can approximate: 
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Therefore, with these approximations, we obtain for the cooling function of HD (Eq. 9) the expression: 
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B H2 molecule 

Concerning the H2 molecules, the analysis is the same than that for the HD molecules with the difference that 
the studied transitions are quadrupolar. Because of the increased number of charged particle regions with molecules 
{H2, HD...), the heating mechanism is different from and more complicated than that of atomic regions. We focus 
our calculations only on the first transition J ~ -^ J — 2. The other transitions are negligible because we consider 
the early phase of gravitational collapse where the matter temperature is below 200 Kelvins which is much lower than 
the transition temperatures for H2 molecules (see Table 1). The Einstein coefficients for the spontaneous emission 
and for the first transition are given by (see Kutner 1984): 
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The relations between the Einstein coefficients are 
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and the radiative density for the transition is: 
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and for the colhsional coefficients in the MaxweU-Boiztmann distribution approximation: 
by anaiogy with the HD moiccuies and with the same ieveis of approximations 
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Thus, from the Eqs (5) and (6), we find for the populations X^'^ and X^ ^: 
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The probability of radiative de-excitation becomes 
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while the molecular cooling for the H2 molecule is given by: 
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C Collision rates 

The colhsional coefficients are given by the expressions: 

Cfo^ = < a.v"° > nH2 (17) 

Cf^ = < CT.i;^^ > n^, (18) 

where a ~ lA^ is the colhsional cross section and nj^^ the density of H2', vhd ^^nd v^^ a-re respectively the velocity 
of HD and H2 given by 



/ 6f£l rn 
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where mn is the hydrogen mass. Notice that, in our model, H2 molecules are more abundant than HD molecules, 
so H2 molecules are the main colhsional species. 



D Total cooling function 

The total cooling function is given by 
where 
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is the ratio between the H2 and HD cooling functions. With the help of the numerical values of the Table 1, we 
deduce: 
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and for the ratio: 
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IV MOLECULAR FUNCTIONS OF A COLLAPSING PROTOCLOUD 

We study a homologous model of spherical collapse of mass M similar to the model adopted in Lahav (1986) and 
Puy & Signore (1996) in which we only consider H2 and HD molecules. 

We have seen (Puy & Signore 1997) from a numerical integration of the cooling functions where twenty rotational 
levels were considered that molecular cooling is important. Moreover we have shown that below 200 Kelvins the 
main cooling agent is the HD molecule. The aim of this paper is to analyse, through our approximations (where the 
only two first excited levels are considered) , the evolution of the cooling and the potentiality of thermal instability. 
First, let us recall the equations governing dynamics of a collapsing protocloud: 
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for the evolution of the matter temperature T„i, 
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for the evolution of the radius r of the collapsing cloud and 
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for the evolution of the matter density n. 

At the beginning of the gravitational collapse the matter temperature increases, then due to the very important 
efficiency of the molecular cooling, the temperature decreases. We consider here this transition regime i.e. the point 
where the temperature curve has a horizontal asymptot (see Puy & Signore 1997): 
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and where the cooling due to molecules can exceed the heating resulting from the collapse. From the equation (22) 
of the evolution of the matter temperature, we deduce 
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Note that the total molecular cooling is negative {hrotai < 0). Moreover, because we consider the early phase of the 
gravitational collapse, we can approximate to a free fall collapse and so the pressure term, 5kTm/2mHr, is negligible 
in the equation (23). Thus, the evolution equation of the radius can be written: 



which leads to 



Moreover, because we know that the evolution of the density is: 
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where Uq and r,, are respectively the density and the radius at the beginning of the collapse. Thus we conclude that 
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We have seen that the regime of temperature that we consider is below 200 Kelvins. Moreover, the abundances 
of HD and H2 remain roughly constant during the early phase of the gravitational collapse (Puy & Signore 1997, 
1998). In this case, the value of the density of HD at the beginning of the collapse is uhd ^ 6 x 10^* cm^''. The 
temperature of the matter is below 200 Kelvins. Thus we can approximate the following terms as follows: 
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These approximations lead to the following simplified expressions for the molecular cooling 
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where A^; Tq and ATq are constants: 

Ao = 2.66 X 10"^^ erg cm"^ K^^/^ s"^ 

To = Tj^o° = 128.6 K 

ATo = t"^ - T^f = 383.4 K 

From our homologous model of gravitational collapse, the density of matter is the density of the species which is the 
most abundant so n = nu2 which gives 
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Finally we find for the total cooling 
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With £_o = 9260, if we take for the mass of the cloud the following approximation: 
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from equations (28) and (32) we deduce for the the radius 
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where the constant 
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FinaUy we obtain for the total molecular cooling 
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V THERMAL INSTABILITY 

We have learned much about thermal instability in the last 30 years since the appearance of the work of Field 
(1965). Many studies have focused on the problem of thermal instabilitie in different situations: formation of galaxies 
(Sofue 1969), expanding medium (Kondo 1970), formation of protostars (Stein et al. 1972), formation of globular 
clusters (Fall & Rees 1985), cooling flows (David et al. 1988) and quasar clouds (Mathews & Doane 1990). 
The correct instability criterion can be derived physically as follows. In an infinite, uniform, static medium of density 
Uo and To, the continuity and force equations are automatically satisfied. If we introduce a perturbation of density 
and temperature such that some thermodynamics variable (pressure, temperature...) is held constant, the entropy 
of the material S will change by an amount 6S, and the heat-loss function, by an amount S'9 = F — A. From the 
equation 
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We know that TdS ~ CpdT in an isobaric perturbation. The corresponding inequality can be written: 
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which characterizes the Field's criterion. The criterion involves constant P because small blobs tend to maintain 
pressure equilibrium with their surroundings when heating and cooling get out of balance. In our case (i.e. at the 
transition regime), a thermal instability could spontaneously be developped if the Field's Criterion is verified: 

dlnA 
dlnTm 
In our case the logaritmic differential of the total molecular cooling function is given by 
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Therefore, the criterion is independant of the mass M of the cloud at the transition regime. In figure 1, we have 
plotted the curve y{Tm) — dlnA/dlnTm- It shows that the Field's criterion is always verified. 
We conclude that a thermal instability is possible at the transition regime, i.e. at dTm/dt = 0, when cooling 
(due to molecules) exceeds heating (due essentially to the collapse). Thus by maintening the same pressure in its 
surroundings, such a blob would get cooler and cooler and denser and denser until it could separate into miniblobs. 
This possibility is very interesting and could give a scenario of formation of primordial clouds. However, a quantitative 
study is necessary to evaluate the order of magnitude of the mass and the size of the clouds. This last point can 
be crucial particularly concerning the composition of the molecular clouds at high redshift such as the Lyman-a 
absorbers in order to evaluate the possible signatures of the primordial clouds. 



VI ON CO MOLECULE 

In the very early Universe the components of the primordial molecular medium arc formed from the elements of 
the primordial nucleosynthesis. Later, as soon as stellar processes occur in proto-objects, other molecular species 
appeared such as CO, CI, HCN . Therefore, their line emission can be seen as an excellent probe of physical 
conditions of very high redshift objects. Among them, the diatomic molecule CO has the strongest lines at radio 
wavelengths, the rotational transitions being spaced by 115 Ghz in the rest frame at z = (this separation, at a 
redshift z, becomes compressed by a factor 1/(1 + z)). Moreover, the carbon monoxide line emission had already 
been detected in four high redshift objects: 

1. C0(3 - 2), C0(4 - 3) and CO(6 - 5) in the ultraluminous IR galaxy F10215+4724 at the redshift z = 2.28 
(Brown and Van den Bout 1992, Solomon et al 1992). 

2. C0(3 - 2), CO{A - 3), C0(5 - 4) and C0(7 - 6) in the cloverleaf quasar at the redshift z = 2.56 (Barvainis 
et al 1994, Barvainis et al 1997). 

3. C0(5 - 4) in the quasar BR1202-0725 at the redshift z = 4.69 (Ohta et al 1996, Omont et al 1996). 

4. C0(5 - 4) in the quasar BRI1335-0418 at the redshift z ^ 4.41 (Guilloteau et al 1997). 

One must also add that, at least for the first three objects, their detectability in molecular transitions may be due 
to high emissivity of the gas and above all magnification from gravitational lensing rather than an extremely large 
mass of gas. Moreover, there is no a priori evidence of gravitational lensing for BRI1335-0418. Therefore, if this 
absence of strong gravitational magnification is confirmed, the fourth object could have a large mass of molecular 
gas which could be associated with the presence of some huge starburst. However, if there is no doubt about the 
reality of the CO detection in these objects, it is more difficult to definitively determine the CO to H2 conversion 
factor in these objects. 

Let us also remark that, here, we study proto structures in the mass range 10^ — 10^° M© i.e. in the mass range 
of dwarf galaxies. Recently, Verier & Hodge (1995) estimated the CO to H2 conversion factor for the very low 
metallicity dwarf irregular galaxy GR8. They found X > 80 x 10^° molecules cm~^ (Km s~^)~^ which is a factor 30 
larger than for our Galaxy {X — N{H2)/Wco where N{H2) is the column density of molecular hydrogen gas and 
Wco is the integration of CO antenna temperature over the velocity profile) . 

Therefore, from these estimations, the calculations of molecular functions could, in principle, be extended to the CO 
molecules for collapsing protoclouds at 2 < 5. The CO molecules could even be a better cooling agent than H2 and 
HD for some temperature ranges. But, in any case -and in particular for the discussion of thermal instability at the 
transition regime- the corresponding total cooling function must be numerically calculated because the excitation of 
many rotational levels must be taken into account. This extended study is beyond the scope of our paper. 
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